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cavity  detachment  point  near  or  at  profile  nose 

cavity  detachment  point  at  profile  trailing  edge 

see  related  nomenclature  below 

see  related  nomenclature  below 

location  of  eigensolution  singularity  on  unit 

circle 

drag  coefficient 
lift  coefficient 

pressure  coefficient  on  wetted  surface 

profile  chord  length,  c  *  1 

element  of  arc  length  in  z-plane 

strength  of  eigensolution 

complex  potential  in  complex  F-plane 

velocity  potential 

stream  function 

stagnation  point  location 

point  at  infinity  in  z-plane 

free-stream  velocity 

intermediate  mapping  complex  plane 

complex  velocity 

intermediate  mapping  complex  plane 
complex  variable  in  the  physical  (x,y)  plane 
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angle  of  atCack 

angular  location  of  stagnation  point  measured 

from  negative  real  axis 

complex  variable  in  unit  circle  plane 

normal  distance  from  profile  chord 

complex  semi -circle  plane,  5  * 

distance  along  profile  chord 

cavity  thickness  at  trailing  edge 

stagnation  point  angular  location  on  unit 

circle  in  5  plane  (cos  y  ■  a/ b) 

angular  location  of  eigensolution  singularity 

on  unit  circle 

value  of  $  at 

value  of  4  at  A2 

0  +  it  complex  logarithmic  hodograph 
flow  inclination 
flow  speed 


Subscripts 

c 


o 


0 


1 


pertains  to  eigensolution  u»c  or  to  the 
cavity  surface 

pertains  to  flat  plate  solution  w0 
and  other  variables  specifically 
associated  with  the  geometric  point  0 
used  on  any  variable  having  zero  as  its 
argument  or  limit 

pertains  to  regular  part  of  solution 
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Introduction 

The  present  paper  bears  upon  the  two-dimensional  inverse  or  design 
problem  for  fully-cavitating  hydrofoils  in  which  one  specifies  the  pressure 
distribution  on  the  profile  wetted  surface  and  then  calculates  that  wetted 
surface  shape  which  will  satisfy  this  prescription.  This  design  problem  is 
certainly  not  new  to  airfoil  designers  and  as  far  as  cavity  flows  are 
concerned,  both  linear  and  nonlinear  design  methods  have  been  worked  out. 

In  the  realm  of  nonlinear  approaches  to  the  present  problem,  the  very 
general  method  of  Yim  and  Higgins  [1]*  is  worthy  of  note  because  it  applies 
to  single  foils  as  well  as  to  cascades  of  profiles  for  all  cavitation 
numbers  in  the  cavity-flow  regime.  Another  approach  has  been  discussed 
superficially  by  Khrabov  [2].  Both  of  these  contain  far  more  generality 
than  is  required  for  this  study  at  zero  cavitation  number.  For  the  direct 
or  off-design  problem  of  exact  cavity-flow  theory,  a  good  example  of  the 
present  level  of  development  is  represented  by  the  work  of  Furuya  [3]  and 
it  is  clear  that  now  one  can  do  both  the  design  and  off-design  problems  for 
fully-cavitating  hydrofoils.  Thus,  one  can  attempt  to  tailor  the  profile 
to  an  entire  set  of  performance  goals  and  failing  that  he  can  at  least 
design  for  the  best  compromise  among  a  set  of  conflicting  requirements. 

According  to  many  authors  [4-7] ,  the  inverse  problem  is  not  thought 
to  present  much  of  a  challenge  at  zero  cavitation  number.  In  this  case, 
the  classical  method  of  Levi  Civita  [7]  can  be  applied  to  an  isolated 
body.  This  view  is  certainly  proper  as  long  as  one  is  content,  after 

^Numbers  in  square  brackets  indicate  citations  in  the  references  listed 
below. 
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prescribing  the  pressure  in  the  circle  plane,  to  accept  whatever  correlation 
between  points  in  the  circle  and  physical  planes  may  result.  Of  course,  such 
a  rudimentary  approach  does  not  lead  to  a  useful  design  procedure. 

The  motivation  for  the  present  investigation  is  that  none  of  the 
literature  on  the  nonlinearized  direct  and  inverse  problems  we  have  surveyed 
so  far  [1-8]  has  made  use  of  the  fact  the  exact  theory  admits  the  existence 
of  a  "point  drag"  or  complementary  function.  While  this  fact  is  of  no 
particular  importance  in  the  direct  problem,  we  have  already  seen  in  the 
case  of  the  linearized  inverse  problem  [9-11]  that  the  complementary  function 
can  play  an  important  role.  For  the  exact  inverse  theory  there  has  been  a 
question  if  a  nonlinear  eigensolution  exists  or  if  it  does,  should  it  be  an 
admissible  component  of  the  solution  [1]?  Therefore,  in  this  study,  we 
explore  these  questions  regarding  the  existence  and  usefulness  of  a  “point 
drag"  or  eigensolution  in  the  nonlinear  theory  under  the  simplest  set  of 
circumstances  and  this  leads  us  naturally  to  the  restrictions  that  the  free 
streamline  flow  pertains  to  an  isolated  profile  and  that  the  flow  be  at 
zero  cavitation  number.  These  simplifications  will  free  us  from  the 
complications  arising  from  non-zero  cavitation  numbers  and  other  boundaries 
in  the  flow  domain  such  as  a  free  surface  or  neighboring  cascade  blades. 

In  this  paper  we  use  the  term  eigensolution  in  the  sense  of  thin 
airfoil  theory  as  suggested  by  the  work  of  Van  Dyke*  because  we  already 
know  that  the  inverse  problem  in  the  theory  of  fully  cavitating  hydrofoils 
is  not  necessarily  unique.  Our  aim  is  to  find  a  sufficiently  weak 

*Perturbation  Methods  in  Fluid  Mechanics,  The  Parabolic  Press,  Stanford,  CA, 
pp.  48-54  (1975). 
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singularity  which  can  be  added  to  the  classical  Levi  Civita  solution  and  can 
then  be  used  to  satisfy  certain  additional  physical  conditions  relating  to 
the  location  of  the  free  streamline  springing  from  the  hydrofoil  nose  and 
thereby  provide  a  unique  inverse  cavity-flow  solution.  After  we  have 
constructed  the  simple  eigensolution,  we  will  examine  some  of  its  properties. 
The  actual  use  of  this  solution  in  the  design  process  will  be  presented 
elsewhere,  although  we  will  start  the  process  so  that  the  potential  usefulness 
of  the  eigensolution  in  design  can  be  seen. 
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Flow  Geometry  and  Conformal  Mappings 

As  noted  above,  this  study  uses  Levi  Civita's  method  [7]  for  the  analysis 
of  the  exact  inverse  problem  for  a  fully  cavitating  hydrofoil  section.  The 
flow  geometry  and  the  principal  quantities  associated  with  the  flow  are 
illustrated  in  Fig.  1.  The  origin  of  coordinates  in  the  z  *  x  +  iy  plane  is 
taken  at  the  stagnation  point  on  the  wetted  surface  of  the  hydrofoil.  This 
point  is  denoted  by  0  in  Fig.  1.  The  chordline  of  the  profile  is  inclined  at 
the  angle  o  with  respect  to  the  x  axis  and  the  free-stream  velocity  U  is  taken 
as  being  parallel  to  this  axis  as  illustrated.  The  flow  separation  point  at 
or  near  the  profile  nose  is  the  point  Aj  as  illustrated  for  a  sharp-nosed 
foil.  The  "upper"  cavity  surface  is  shown  as  the  dashed  curve  extending  from 
Aj  to  the  point  o'  at  downstream  infinity.  In  the  case  of  a  round-nose 
profile,  Ai  can  lie  on  the  upper  wetted  surface  behind  the  leading  edge.  This 
case  is  not  illustrated  in  Fig.  1.  The  point  A2  denotes  the  location  of  the 
trailing  edge  of  the  wetted  surface.  The  lower  surface  of  the  cavity  leaves 
the  wetted  surface  at  A2  and  extends  as  shown  by  the  dashed  line  to  the  point 
o’  at  downstream  infinity. 
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Let  the  coordinates  of  a  typical  point  on  the  wetted  surface  be  denoted 

by  x  and  y  and  those  on  the  upper  cavity  by  x^  and  y^  as  shown  in  Fig.  1. 

While  the  orientation  of  the  profile  in  the  z  plane  is  convenient  for  the 

purposes  of  analysis,  the  x-y  system  is  not  always  a  convenient  reference 

frame  for  foil  and  cavity  contours.  For  this  purpose  we  use  a  coordinate 

system  with  the  abscissa  along  the  chordline  as  shown  by  the  distance  a 

measured  from  the  profile  nose.  The  ordinates  of  the  wetted  surface  are 

then  given  in  terms  of  a  as  n(o)  and  the  upper  cavity  ordinates  are  given 

by  n  (a).  At  the  trailing  edge  of  the  profile,  the  cavity  thickness  is 
c 

n  *  T.  These  quantities  are  also  shown  in  Fig.  1.  In  the  a,n  system 
c 

the  stagnation  point  0  is  located  at  (°0»h0)  as  illustrated.  The  trans¬ 
formation  between  the  (x,y)  and  (c,n)  systems  is 

o  +  in  *  a  +  in  +  ze  ,  (1) 

o  o 

where  z  is  the  complex  variable,  z  =>  x  +  iy,  and  a  is  the  angle-of -attack 
as  measured  by  the  inclination  of  the  chordline  with  respect  to  the  x  axis 
and  free-stream  velocity  U. 

The  conformal  mappings  start  with  the  complex  potential  in  the 
z -plane , 

F  =  <j)  +  i\j»  ,  (2) 

whc.-**  d>  is  the  velocity  potential  and  \p  is  the  stream  function.  As  is 
customary,  we  adjust  these  quantities  to  make  $  =  0  at  the  stagnation 
point,  0.  The  stream  function  is  taken  to  be  zero  all  along  the 
stagnation  streamline.  Therefore,  the  boundaries  of  the  flow  can  be 
represented  by  a  cut  all  along  the  real  axis  in  the  F-plane  as  shown  in 
Fig.  2.  Note  that  the  wetted  surface  extends  from  the  stagnation  point 
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at  0  to  the  trailing  edge  at  A2  and  the  lower  cavity  surface  from  A2  to  O' 
must  lie  along  the  lower  surface  of  the  cut.  This  is  so  because  downstream 
from  the  stagnation  point,  the  velocity  potential  increases  in  the  flow 
direction  and  the  stream  function  decreases  outwardly  from  the  foil  or 
cavity  surface.  On  the  arc  QAi,  the  flow  direction  is  reversed  with  a 
consequent  reversal  in  the  gradients  of  $  and  1J1  so  that  the  point  A}  is 
on  the  upper  edge  of  the  cut. 

One  can  use  the  mapping, 

W  -  ✓?  (3) 

in  crder  to  map  the  F-plane  outside  the  cut  into  the  upper  half  of  the 
W-plane.  Corresponding  points  are  shown  by  the  locations  of  o',  Aj  and 
A2  and  o'  in  Fig.  3.  As  before,  the  cavity  surfaces  are  shown  as  dashed 
lines  and  the  wetted  surface  is  shown  by  the  solid  line.  Let  the  values 
of  $  at  Aj  and  A2  in  the  F-plane  be  $1  and  $2  respectively.  Then  these 
points  are  at  W  *  and  W  ■  -  end  the  midpoint  of  the  distance 
between  A}  and  A2  is  located  at 


The  distance  between  this  midpoint  and  Aj  is 


/?7  + 

b - - - -  .  (5) 


These  distances  are  also  shown  in  Fig.  3. 
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One  can  now  use  the  mapping 


W  - 


(Z  +  cos  y)  , 
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(6) 


where 


► 

I 


a  '*1  '  '*2 
COS  Y  -  yT  ■  - - — 

+ 


(7) 


in  order  to  map  Aj  into  Z  *  +  1  and  A2  into  Z  ■  -  1.  The  point  0  maps  into 
Z  *  cos  y  as  shown  in  Fig.  4.  The  significance  of  the  point  C,  also  shown 
in  this  figure,  is  discussed  later.  Then  the  Joukowski  transformation, 


z  -  -  7  U  +  . 


(8) 


maps  the  upper  half  of  the  Z-plane  onto  the  interior  of  the  unit  circle  in 
the  ? -plane.  The  inverse  of  this  mapping  mist  be 

C  -  -  Z  + V  Z2  -  1  (9) 

in  order  to  make  the  point  o'  map  into  the  origin  of  the  c -plane.  Various 
corresponding  points  are  denoted  by  Aj,  A2  and  0  in  Fig.  5.  The  arc  of 
the  unit  semi-circle  corresponds  to  the  wetted  surface.  Coordinates  on 
the  wetted  surface  are  given  by  5  *  e^-S  as  shown.  The  stagnation  point  is 
at  B  ■  y*  The  upper  surface  of  the  cavity  is  on  the  real  axis  between 
A^  and  0*  and  the  lower  surface  of  the  cavity  is  on  the  real  axis  between 
A2  and  o'  as  marked  in  Fig.  5.  One  can  now  use  Eqs.  (3)  through  (8)  in 
order  to  write  the  composite  of  the  preceding  mappings  as 


&&&& :  ■  *  ; 
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F  -  b2[cos  y  -  y  (t  +  “)]' 


Then  we  Introduce  the  complex  velocity, 


w  -  u  -  iv  -  -7—  , 

dz  ’ 


in  order  to  write 


-i0  .  2 r 

w  -  qe  *  b  [cos  y 


1  {  1%  1  d  ? 

-2  +7)](-5  - 


wdz  -  dF  -  b2  4  (t  +  7)  -  cos  r]  (t  -  7)  T5-  • 


These  quantities  are  now  used  to  define  the  logarithmic  hodograph  or  w-plane 


-iai(c)  \  dF  w  _  £  -i0 

5  U  dz  U  U 


exp[in  -  i0] 


Therefore  we  have 


«(t)  ■  8  +  iin  ^  »  0  +  ix  , 


where  x  ■  in  j*-.  On  the  free  streamlines  q  ■  U  so  that  x  *  0  there.  In  the 
C -plane,  these  free  streamlines  are  on  the  real  axis  and  at  o'  we  know  that 


0  *  0  also.  Therefore, 


u(0)  -  0 


and  u(c)  is  real  when  ;  is  real.  At  the  stagnation  point  q  +  0  so  that 
x  ♦  -  •  there.  The  flow  directions  differ  by  ir  on  either  side  of  0  and 
so  the  u»-plane  with  the  various  corresponding  boundaries  can  be 
represented  as  illustrated  schematically  in  Fig.  6. 
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We  can  new  use  Eqs.  (13),  (14)  and  (15)  in  order  to  write 

.2  -T+i0  ,  , 

dz  «  dx  +  idy  -  — — -  [j  (£+•£■)-  cos  y](;  -  -£■)  .  (17) 

IB 

On  the  wetted  surface  £  *  e  and  Eq.  (17)  leads  to 

2b2  -t 

dx  *  — jj—  e  [cos  y  -  cos  g]  sin  g  cos  0  dg 

and  -  .  (18) 

2b2  -T 

dy  *  e  [cos  y  “  cos  g]  sin  g  sin  0  dg 

Note  that  dy/dx  ■  tan  0  as  it  should  if  the  wetted  surface  is  to  be  a  stream¬ 
line.  On  the  upper  surface  of  the  cavity  t  ■  0  and  arg  £  -  it  so  that  Eq.  (17) 
leads  to 

2 

dxc  "  ¥"  cos  8  +  £)  "  cos  y](c  - f5- 

(19) 

dyc  -  §-  sin  0  [|  (£  +  £)  -  cos  y](c  ~ 
provided  that  -  1  <  (  <  0. 

Returning  to  Eq.  (17),  we  can  write  the  square  of  the  arc  length  along 
the  wetted  surface  as 

2  —  h2  — r  2 

(ds)  «  dz  dz  *  {2  y-*  e  [cos  y  -  cos  g]  sin  g  dg}  , 


where  we  have  also  made  use  of  Eq.  (18)  because  it  applies  to  the  wetted 
surface.  But  from  Eq.  (15)  we  have 
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T**tn^-»ln/l-C  , 

U  P 

so  chat 

.2  _ _ 

ds  -  2  [(cos  y  -  cos  8)  sin  0//1  -  Cp(8)JdB  (20) 

on  the  wetted  surface  of  the  profile.  As  we  have  noted  previously,  the  flow 
directions  differ  by  s  on  either  side  of  0.  Therefore,  if  the  sign  of  dz  is 
positive  on  the  arc  OA2,  it  will  be  negative  along  arc  0A}.  As  a  result  of 
this  difference  ds  might  have  a  like  sign  change  in  these  two  regions.  Just 
how  this  might  occur  depends  on  the  form  of  /l  -  Cp(8)  in  any  particular 
case.  Therefore,  we  will  defer  consideration  of  this  question  to  a  later 
place  in  the  development. 

Hydrodynamic  Forces 

The  development  of  general  formulae  for  the  hydrodynamic  forces  on  the 
profile  depend  upon  certain  properties  of  the  function  u(^)  which  result 
from  the  previously  noted  fact  that  «({)  is  real  when  c  is  real.  For  then 
one  can  apply  Schwartz's  principle  of  symmetry  in  order  to  write  w(c)  -  u»( <;) 
and  thereby  obtain  the  analytic  continuation  of  u(c)  into  the  lower  half  of 
the  unit  circle  [6] .  Thus  we  can  write  for  a  prescribed  modulus, 

r0  "  kl  <  1» 

9(8)  -  ix(S)  -  9(-6)  +  ix(-8)  , 

or 

T(-8)  -  -  t(8) 

and 


6(“8)  -  9(8) 
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Hence  t  is  an  odd  funcClon  of  f$,  or  Im  (,  and  6  Is  an  even  function  of  8,  or 
Im  (. 

Using  the  result, 

eT  -  /I  -  C  (21) 

P 

as  noted  above,  and  the  fact  that  w(()  is  now  defined  inside  the  unit  circle 
one  finds  [7]  from  the  calculus  of  residues  that 

.2 

CL  "  I  C4“,<0)  cos  Y  -  u>,,(0)]  (22) 

and 

S*T4t“'<0)]2  •  <23> 

where  the  quantity  c  is  the  profile  chord  which  is  taken  as  unity  in  this 
work.  The  moment  can  be  calculated  after  the  complete  solution  has  been 
found. 

The  Form  of  oi(c)  Near  the  Stagnation  Point 

This  is  also  a  well  known  result  which  we  shall  review  briefly.  The 
form  of  u>  near  0  is  dominated  by  the  fact  that  on  a  smooth  contour  u  jumps 
by  the  amount  ir.  In  particular,  as  one  traces  the  profile  surface,  starting 
at  A}  in  Fig.  1  and  then  passes  through  0  while  proceeding  to  A£,  the  jump 
is  a  decrease  in  0.  This  is  precisely  the  behavior  exhibited  by  the  real 
part  of  the  analytic  function  ifcn((  -  e  T).  However  this  function  by  itself 
does  not  have  the  proper  symmetry  needed  for  admissible  forms  of  at.  But  if 
we  subtract  from  it  a  similar  function  which  has  a  like  jump  at  the  image  of 
0  with  respect  to  the  real  axis  in  the  (-plane,  we  preserve  the  necessary 
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behavior  ae  0  and  also  satisfy  the  preceding  symmetry  requirements.  Aside 
from  arbitrary  additive  constants,  this  function  has  the  form 


Finally,  we  require  the  ia( 0)  -  0.  Because  of  this  condition  the  resulting 
function  which  provides  the  flat  plate  solution  is  [4] 


u  (c)  -  0  +  it 

o  *  o  o 


-  o  +  iin 


1  -  ?eJ 


with  y  »  ir  -  a  in  the  case  of  an  isolated  flat-plate  profile.  For  this  case, 

one  can  show  that  this  function  has  the  flow  direction,  6  »  -  a  behind  the 

’  o 

i  6 

stagnation  point  or  9q  ■  ir  -  a  ahead  of  the  stagnation  point  on  c  ■  e  . 
Therefore,  the  wetted  surface  in  the  z-plane  is  a  straight  line  through  0 
with  its  trailing  edge  Inclined  at  the  angle  -  a  with  respect  to  the  positive 
real  axis.  Moreover,  since  t0  vanishes  on  the  real  axis  in  the  c-plane,  the 
free  streamlines  have  -  0  as  required.  Thus  we  can  write  Eq.  (24)  as 

1y 

»0<C)  -  *H(B  -  y)  -  a  +  Un  »  (24a) 

where 


H(8  -  y) 


o  ,  e  <  y 

1  ,  3  >  Y 


is  the  Heaviside  function 
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The  contribution  of  u  to  C.  and  C_  follows  from  Eqs.  (22)  and  (23). 

o  L  D 

om  Eqs.  (22),  (23)  and  (24)  we  get 


■  2it  —  sin  a  cos  a 
L  Uc 

o 


C  *  2it  —  sin  a  (26) 

D  Uc 

o 

and  we  see  that  L/D  *  cot  a  as  is  proper  for  the  flat  plate.  We  can  also  use 

2  T 

the  relationship  C  ■  1  -  e  to  find  the  pressure  distribution  on  the  plate. 
The  result  is 


4  sin  8  sin  y 


po  (cos  8  -  cos  y)  +  (sin  6  +  sin  y)‘ 


From  this  result,  we  see  that  when  8  *  Y»  CD  ■  1  and  when  8  “  0  or  i, 

*o 

Cp  »  0  as  required. 

*o 

Continuing  the  study  of  the  flat-plate  solution,  we  can  rewrite  Eq.  (24) 
for  points  on  the  wetted  surface  as 


u>  *  ir  -  a 
o 


+  Y  +  iin[(eiY  -  ei0)/(e  lY  -  eiS)]  , 


in  which  case 


ir  _  ie  _ 

- - - —  -  /i  -  c 

e’iY  -  eiS  Po 


After  some  manipulation  we  find  that 


✓  l  -  C 


2  cos  g  -  cos 

_  „\ 2  ,  r 


Po  (cos  y  ~  cos  g)  +  (sin  y  +  sin  g) 

Next  we  can  put  g  ■  w  -  5  and  y  =  it  -  o  in  this  expression  and  in  Eq.  (17) 
when  it  is  written  for  points  on  the  wetted  surface.  That  is, 

2 

,  b  (cos  5  -  cos  a)  iQ  r/-  _  -*2  .  r 

dz  -  -  7= — i - ? - r  e  l  cos  5  -  cos  al  +  (sin  5  +  sin  a]  sin  5  d£ 

U  | cos  §  -  cos  a|lA  1  J 

2 

A8  b  2 

=*  -  e  sgn(cos  £•  “  cos  a)  2  —  [sin  £  +  sin  a  sin  £  -  cos  a  cos  £  sin  £]d£ 

19 

We  can  dispose  of  the  product  -  e  sgn  (cos  £  ~  cos  o)  by  observing  that  on 

the  wetted  surface  when  0  <  £  <  a  the  flow  direction  is  0  *  it  -  a  and  when 

1 9  “1  oi 

a<E<x,0*-a.  Therefore  in  the  first  of  these  cases,  e  -  -  e  and 
1 9  “1  ct 

in  the  second  e  *  +  e  .  Hence  in  either  of  them  the  product 

-ei9  sgn  (cos  5  -  cos  a)  ■  e  ia.  Next  we  can  introduce  d(o  +  in)  *  eiadz 

from  Eq.  (1)  with  the  result  that 


b  2 

d(a  +  in)  ■  2  —  [sin  £  +  sin  o  sin  £  -  cos  a  cos  £  sin  £]d£ 

This  last  result  implies  that  n  ■  0  as  is  proper  for  a  flat  plate  and  do  is 
simply  the  arc  length  ds  along  the  wetted  surface  measured  from  the  profile 
nose  where  a  *  0.  We  can  integrate  this  last  equation  from  o  *  0  (£  *  0) 
to  some  value,  0<o<l(0<£<ir),  and  get 


8  "  2  IT-  [(  1  -  cos  ?)(1  * 


[1  +  cos  5])  +  (5  -  sin2  5)  -i°-1  . 
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The  profile  is  to  have  unit  chord  however,  so  that  when  5  «  ir,  o  *  s  *  1 
and 


U  4  +  ir  sin  o  v 

When  Eq.  (29)  is  used  in  Eqs.  (25)  and  (26),  one  obtains  the  well  known 
Rayleigh  formulae  [12]. 

Additional  properties  of  the  flat-plate  solution  which  are  important 

for  the  present  considerations  relate  to  the  shape  of  the  upper  surface  of 

the  cavity.  In  Fig.  5,  £  will  be  on  the  negative  real  axis  for  points  on 

this  part  of  the  cavity.  Therefore,  let  us  put  C  *  -  C  where  £  is  a  real 

c  c 

i  Cl 

positive  number.  Then  we  can  use  d(o  +  in)  *  e  dz ,  and  write  Eq.  (17)  in 
the  form 


d(o  +  in)  -  jp  exp  i[a  +  0oUc)]d{[j  (c  +•£“)-  cos  a]2}  . 


But  now 


exp  i[a  +  ej  -  e'ia(eia  -  Cc)/(e_ia  -  £c) 

and  so  we  have 

2 

d(o  +  in)  =  j  [Cce‘ia  -  2  +  (21  sin  a)/?c  +  2/^  -  eia/c2]dcc  . 

The  integration  in  this  case  starts  at  A,  where  C  *  1  and  a  -  n  *  0  and 
proceeds  to  some  value  of  oc>nc  corresponding  to  1  <  ?  <  0.  This  leads  to 
the  parametric  representation  of  hc(oc)  in  terms  of  which  is  given  by 
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and 


o 

c 


ibiiilicl! 

2  U  Cc 


1  •  ?c 

- j-S-  +  2  ine 

-  K 


(30) 


(31) 


Equation  (29)  gives  the  value  of  jj—  to  be  used  in  Eqs.  (30)  and  (31). 

It  is  useful  to  find  the  ordinate  of  the  upper  surface  of  the  cavity 
above  the  trailing  edge  of  the  wetted  surface.  This  can  be  done  by 
calculating  the  quantity 


(i-?cr 


(i  -ccr 

2c 


cos  a-2 


for  several  values  of  Cc  and  prescribed  values  of  a.  Then  when  =  1,  one 

can  plot  contours  of  a  *  constant  in  an  f  -  plane  and  note  for  each  value 

of  a  the  value  of  c  corresponding  to  f  *  2(4  +  it  sin  a).  These  values  of 

c 

C  can  then  be  used  in  Eq.  (31)  in  order  to  compute  the  value  of  n  (1)  for 
c  c 

each  value  of  a.  It  was  found  when  these  points  were  plotted  in  an  nc(l)  - 
plane  that  a  linear  relationship  fits  the  data  for  0  <  a  <  10°.  When  a  is 
measured  in  degrees,  this  line  has  the  equation 

n  (1)  =*  .0294a°  .  (32a 

c 

The  corresponding  relationship  for  a  measured  in  radians  is 

n  (1)  -  1.684a  .  (32b 

c 

The  computed  results  are  tabulated  below. 
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Table  1.  Cavity  Thickness  at  the  Trailing  Edge  of  a  Flat  Plate 


a0 

Cc(l) 

nc(0 

0.0 

-  _ 

.0000 

0.5 

.1712 

.0147 

1.0 

.1707 

.0297 

2.0 

.1698 

.0585 

3.0 

.1689 

.0879 

4.0 

.1680 

.1174 

5.0 

.1669 

.1463 

7.0 

.1650 

.2049 

10.1 

.1620 

.2939 

These  results  will  be  used  in  order  to  start  the  functional  iteration  after 
the  formulation  of  the  theory  has  been  completed.  As  a  final  note,  we 
observe  that  the  cavity  thickness  at  the  trailing  edge  of  a  flat  plate 
according  to  linearized  theory  [13]  is 


nc(l)  «  1.681a  , 


which  can  be  compared  to  the  corresponding  expression  from  Eq.  (32b).  Within 
the  range  of  attack  angles  considered  here  the  trailing-edge  cavity  thickness 
at  zero  cavitation  number  is  about  the  same  when  estimated  by  linearized 
or  nonlinear  theory.  For  larger  values  of  attack  angle,  we  would  expect 

estimates  from  linearized  theory  to  exceed  nonlinear  cavity  thicknesses. 

The  flat-plate  function  u>0(5)  is  traditionally  considered  to  possess 
all  of  the  singular  behavior  of  the  function  u>(c).  The  shape  of  the  smooth 
body  is  then  represented  by  an  analytic  function  011(5)  which  is  regular 
inside  and  on  the  unit  circle.  It  must  also  satisfy  the  same  symmetry 
requirements  as  are  imposed  upon  u)q  and  we  must  also  insist  that  oj^(0)  *  0. 

Then  one  traditionally  puts  a >(?)  ■  u)q  +  As  mentioned  previously,  we 

will  add  to  this  customary  sum  a  new  function,  0^(5),  which  is  the  analog 
of  the  point-drag  function  of  linearized  theory.  We  will  now  explore  the 
properties  of  this  eigensolution. 
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A  Simple  Eigensolution 

The  complementary  function  ^  “  ®c  +  *TC»  *s  t0  ^  determined  from  the 

requirements  that  T  *  0  on  the  cavity  and  the  foil,  8  =0  along  the 

c  c 

stagnation  streamline  and  that  o>  vanishes  at  infinity.  A  function  of  the 

c 

form  which  satisfies  these  conditions  can  be  found  most  easily  by  considering 
the  flow  in  the  F-plane,  Fig.  2.  If  we  take 


01 

c 


0  +  it  -  — 

c  C 


(33) 


where  E  is  a  real  constant,  we  have  a  function  which  satisfies  the  necessary 
requirements.  The  two  conditions  *  0  and  q^  =  0  on  both  the  cavity  and 

the  profile  wetted  surface  can  be  satisfied  by  any  member  of  the  family  of 
functions  having  the  form  F-™,  0  <  m  <  1.  But  the  condition  0^  =  0  along  the 
entire  stagnation  streamline  can  be  satisfied  only  when  m  *  1/2.  This  choice 
for  the  complementary  function  seems  to  offer  the  advantage  that  it  will 
cause  less  alteration  of  the  upstream  flow  field  inclination  than  other 
possibilities.  Moreover,  it  is  the  only  choice  from  amongst  the  functions 

— ni 

F  which  gives  the  correct  branching  of  the  flow  and  it  appears  to  be  the 
most  convenient  choice  for  further  analysis.  Consequently,  we  shall  adopt 
this  functional  form  for  the  simple  eigensolution  in  this  work.  The  word 
simple  indicates  that  the  branch  point  for  this  solution  is  coincident 
with  the  stagnation  point  at  0  in  Fig.  1.  We  shall  generalize  this  result 


later 
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It  has  been  noted  by  0.  Furuya*  that  Eq.  (33)  is  just  the  single-spiral- 
vortex  function  proposed  by  Tulin  as  a  useful  representation  for  cavity 
termination  in  the  direct  problem  at  non-zero  cavitation  numbers  [14].  The 
small-scale  structure  of  this  function  is  responsible  for  its  name  as 
discussed  by  Tulin  on  p.  21  of  Ref.  [14],  In  connection  with  the  present 
application,  we  may  also  mention  Tulin's  double  spiral  vortex  model  [14] 
which  we  shall  write  as 

to  =6  +  it  =  D  in  /F  =  D[ln  /r  +  i(*r&  ■?■)]  , 

cd  Cd  Cd  1 

where  D  is  a  real  constant  and  r  =  | F | .  At  points  in  the  physical  plane 

which  are  far  removed  from  the  profile  and  the  cavity,  F  ~  U  and  ui  is 

Z  °d 

not  bounded  at  the  point  0*.  As  it  stands,  this  form  of  the  double-spiral 

vortex  violates  Eq.  (16)  and  it  will  generally  not  produce  a  null  pressure 

coefficient  everywhere  on  the  wetted  surface.  Therefore,  it  is  not  an 

admissible  candidate  for  an  eigenfunction.  This  is  not  to  say  that  other 

logarithmic  forms  for  an  eigenfunction  can  not  be  acceptable.  This  writer 

has  not  found  one  as  yet,  however.  Therefore,  we  shall  content  ourselves 

with  the  form  of  prescribed  by  Eq.  (33)  and  restrict  the  present  analysis 

to  eigensolutions  of  this  form. 

Equation  (10)  can  be  used  to  represent  u>  in  the  c -plane  as 

c 


*Private  Communication  (March  1985). 
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In  the  c -plane  a>c(0)  -  0  and  when  ?  is  real,  (j»c  is  real.  Moreover,  <uc(c) 

is  an  analytic  function  which  is  regular  inside  the  unit  circle  and  which 

— iy 

has  simple  poles  at  {  «  e  .  Note  from  Eq.  (34)  that  on  the  "nose  cavity 

0c  0+  as  c  -*■  -  0  and  on  the  "tail  cavity"  0  ♦  0-  as  5  ■#>  +  0,  as 

-io> 

illustrated  in  Fig.  1  also.  From  Eq.  (14),  w^  ■  Ue  c  and  we  see  that 

the  structure  of  oi  leads  to  an  essential  singularity  in  w  at  the 
c  c 

stagnation  point  0.  The  complex  velocity  w^  is  bounded  at  this  point 
however*,  and  a  smooth  foil  contour  will  pass  through  z  ■  0  as  will  be 

ig 

shown  below.  At  points  on  the  unit  circle,  5  -  e  ,  we  have 

t  -  0  (3 

c 

and 

6c(8)  -  b (cos  „  -  cos  el  ’  2b  sln  ZTJ  sln  •  <3 


*For  points  on  the  unit  circle  very  near  the  stagnation  point,  8  *  y»  let 
0  ■  y  -  e  with  e  «  1.  Then 

•  E 

<*>c<Y,e)  - - - - 

eb[sin  Y  +  y  cos  y] 

and  to  0(e  *)  w£  -  exp  -  iE/eb  sin  y*  Therefore  | w^ |  «  q^  <  1  as  e  +  0. 
Next  consider  an  interior  point,  ?  -  re1^.  Now  let  r  -  1  -  p  with 
0  <  p  «  1.  Then  to  0(p  ),  u>c(y,p)  ■  — g-  ain  t  ■  iinqc/U.  Consequently, 
qc/U  ■  exp  -  E/pb  sin  y  and  qc/U  +  0  as  p  +  0.  Appropriate  linear 
combinations  of  these  two  cases  can  be  used  to  consider  other  limiting 
paths  but  these  give  no  new  information  about  the  boundedness  of  q  /U. 
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Since  t  ■  0  on  the  wetted  surface,  we  expect  u  to  make  no  contribution 
c  c 

to  the  lift  although  the  singularity  at  the  stagnation  point  should  lead  to 
a  drag  force.  Making  use  of  Eqs.  (22),  (23)  and  (34),  we  find  that 


-  0 


and 


From  Eq.  (36)  we  see  that  the  flow  direction  is  not  defined  at  the 
stagnation  point,  8  =»  y.  If  g  <  y  however, 

0  -  -  —r - - r  <  0 

c  b[cos  8  -  cos  y] 

along  the  arc  OA2.  If  g  >  y 

e  .  — - - - —  >  0  , 

c  b[cos  y  -  cos  g] 

along  the  arc  OA^.  Therefore  8^  changes  sign  at  g  *  y  and  as  noted 
previously,  the  real  part  of  u 1  has  a  jump  of  x  as  it  passes  through  0. 

In  view  of  Eq.  (15)  which  requires  that  x  -  Znq/U,  this  situation  suggests 
a  strong  similarity  between  Fabula's  "step-profile”  solution  for  the 
linearized  theory  [15]  and  the  present  eigensolution.  The  present 
application  and  those  of  Refs.  [13]  and  [15]  are  different,  however. 


Simple  Eigensolution  Geometr 
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If  we  put  t  ■  E//F,  we  find  for  a  profile  of  unit  chord  that 


2E  1  it  , 
—  dt 


^Vltd' 

'  t3 


which  has  the  indefinite  integral. 


2  it  it  it 

E  e  .  .  e  .  f  e _ 

z  ’  —  rr ♦ 1  t  + !  ~r dt  • 


Completion  of  this  integration  can  be  carried  out  in  four  parts  starting  from 
either  side  of  the  stagnation  point  where  z  *  0  and  t  +  *. 

For  example,  on  the  arc  OAj  we  recall  that  8  >  Put  t  ■  tj,  Y  ■  »  ~  5 


and  »  tr  -  5l  then 


m  a _ 

1  "  b(cos  -  cos  6) 


Because  0  <  <  6,  we  have 


b(l  -  cos  6)  *  C1  * 


After  separating  Eq.  (37)  into  real  and  imaginary  parts,  we  find  for  this 
part  of  the  wetted  surface,  the  coordinates 


2  cos  t,  sin  t, 

-r  —r1  -  -tt1  +  “(■,) 


2  sin  t,  cos  t 


__1  +  __ — L  +  Si(t  )  -1 
2  t,  ^  1'  2 
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On  the  rest  of  the  wetted  surface,  that  is,  on  arc  OA2  between  the  stagnation 
point  and  the  tail,  we  know  that  8  <  y  and  we  put  t  -  -  t£,  y  ■  ir  -  6  and 
02  ■  *  "  §2*  Then 


t 


2 


_ E _ 

b[cos  6  -  cos  £2] 


In  this  case,  6  <  £2  <  ir  and  it  follows  that 


<  < 


2  b(l  +  cos  6) 


Then  we  have 


and 


X2  “  U 


cos  t„  sin  t„ 

— r2  -  -E72  ♦ 


sin  t.  cos  t. 

2  +  +  Sl(t2, . , 


(39) 


On  the  cavity  surface  A^O  the  integration  starts  at  £  -  -  1  and  ends 
somewhere  between  £  »  -  1  and  ?  »  0.  Note  that  the  value  of  t  at  {  »  -  1 
equals  that  of  tj  when  8*0.  So  that  if  we  put  t  *  t3,  y  *  ir  -  6  and 
£  *  -  £3,  we  have 

E 


b[{  (?3  +  j-)  - cos  y] 


then  the  upper  cavity  joins  the  nose  arc  of  the  wetted  surface  and  its 


coordinates  are 
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_2  cos  t,  sin  t, 

*3-1-  -sp - E72  ♦  “('j) 


g2  sin  t3  cos  t3 


In  this  case 


0  <  t  <  -TT-. - Tv*  . 

c  b(  1  -  cos  6) 


Finally,  the  cavity  surface  from  the  trailing  edge  is  obtained  by 


putting  t  *  -  t4,  y  ■  it  -  6  and  c  *  C 4,  where  0  <  54  <  1.  Therefore 


b[cos  6  - 


0  <  C4  <  b(l  +  cos  6)  ' 


Then  we  have 


2  cos  t,  sin  t, 

*4-!- 

C4  4 


_2  sin  t. 
E  U 


cos  t , 


+  Si[tJ  -  7 


As  we  found  for  the  upper  cavity,  the  lower  cavity  surface  starts  smoothly 


from  the  wetted  surface. 


■  A V jV V*V>s--> V-1-* 
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Equations  (38),  (39),  (40)  and  (41)  will  provide  the  shape  of  the 
wetted  surface  and  the  cavity  surfaces  for  the  simple  complementary  function 
of  Eq.  (34).  These  equations  contain  the  undetermined  ratio  E^/U,  however. 
We  shall  consider  the  ratio  E/b,  which  determines  the  strength  of  the 
complementary  solution,  as  a  parameter  which  we  can  prescribe  —  at  least 
for  the  time  being.  We  also  consider  the  value  of  y  (or  6)  to  be  known. 
Therefore,  we  need  to  “scale"  our  results  in  order  to  obtain  a  profile  of 
unit  chord.  Since  we  anticipate  that  the  complementary  function  can 
produce  a  rounding  of  the  wetted  surface  nose,  the  scaling  procedure  must 
account  for  this  possibility.  Accordingly,  we  shall  need  to  determine 
explicitly  the  location  of  the  apex  of  the  wetted  surface  nose  with  respect 
to  the  profile  chordline. 

Figure  1  shows  the  chordline  for  the  sharp-nose  or  for  the  round-nose 
case  when  the  upper  cavity  separates  at  the  leading  edge  with  respect  to 
the  hydrofoil  chord.  The  geometry  for  the  rounded  nose  with  the  separation 
point  on  the  upper  wetted  surface  behind  the  apex  of  the  wetted  surface 
contour  is  Illustrated  in  Fig.  7.  The  apex  is  located  at  the  origin  of 
o,  n  coordinates  in  this  illustration  with  the  n  axis  being  tangent  to 
the  contour  at  this  point.  Denote  the  x,  y  coordinates  of  this  point  by 
xa»  Ya»  Then  since  the  o-axis  is  normal  to  the  n  axis,  we  see  that  at  the 
apex  the  slope  of  the  contour  is 

dy 

■  tan  (y  -  a)  -  cot  a  .  (42) 

a 

We  will  restrict  our  attention  to  those  cases  in  which  the  apex  is  on  the 
arc  OAj .  Let  t  -  tj  -  ta  at  the  apex.  Then  from  the  equations  preceding 
Eq.  (38),  which  define  tj,  we  have 
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Ca  “  b(co8  5a  -  cos  5)  *  0  *  5a  <  6  »  (43) 

and  Eq.  (43)  can  Chen  be  used  In  Eqs.  (38)  Co  define  xa  and  ya  once  Che  value 
of  ?a  (or  8a)  has  been  found.  Thus,  we  muse  deCermine  Che  unknowns  E^/U  and 
?a  in  Cerms  of  Che  prescribed  quancicies  E/b  and  6(or  y).  Two  condiCions 
are  available  for  Chis  purpose.  The  firsc  is  given  by  Eq.  (42).  The  second 
will  be  chac  'he  profile  has  a  unic  chord. 

An  alCemaCe  form  of  Eq.  (42)  is 

cos  a  *  sin  c 


sin  a  ■  cos  c 


which  follows  from  Che  complex  equacion  jusC  above  Eq.  (37)  when  we  puC 
c  *  Ca.  We  will  use  Eq.  (43)  as  Che  appropriaCe  expression  for  Che  slope 
of  Che  foil  concour  aC  Che  apex.  Now  leC  us  dif ferenciace  Eq.  (1)  so 


d(o  +  in)  -  e  dz  . 


Then  from  Che  complex  equacion  jusC  preceding  Eq.  (37)  we  have 


e  *ad(®  +  in) 


2E  1  ic.„ 

—73 e  dc 


Starring  from  0  where  (o,n)  “  (o0,n0)»  **  incegrace  Eq.  (45)  Co  A2,  where 
(o,n)  *  (1,0).  This  seep  gives 
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Eliminating  the  sum  oQ  +  iq0  from  these  two  equations,  we  get 


j  cos  a  -  F  -  f(tj 


U  r  \ 

—  sin  a  -  G  -  g(t  ) 


In  Eqs.  (46)  we  have 


x  U  cos  t  sin  t 

F  -4 - 


y  U  sin  t  cos  t 

G  *  - T"2-  +  ~t~ *  Si(tt) 


which  contain  known  quantities  because  tt  is  known.  The  remaining  pair  of 
functions. 


cos  t  sin  t 

f(ta)  *-VL--r-t  +  clltJ 

t  a 

a 


g{ta) 


sin  t  cos  t 
a  ,  a 


r*- — r-* ♦  “(«.)  -i 


contains  ta  which  depends  upon  the  unknown,  (Ja  *  n  -  5a*  Thus,  Eqs.  (46)  are 
two  simultaneous  equations  containing  the  quantities  U/E2,  a  and  6a  which 
must  be  determined.  Therefore,  Eqs.  (44)  and  (46)  form  a  determinate  system 


kThe  F  introduced  here  is  not  to  be  confused  with  F  *  $  +  iiji  from  Eq.  (2). 


A  v.  A  A 
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which  can  be  solved  by  iteration.  In  order  to  do  this,  we  can  write  the 
complete  system  as 


G  -  «(C  ) 
ta"  °  f(t-y 

a 


(49) 


cot  t  *  tan  a  (50) 

a 

and 

F  -  f(t  )  G  -  g(t  ) 

U  - - L_  - - - (51) 

jy2  cos  a  sin  a 

In  the  derivation  of  this  system  we  have  assumed  that  the  apex,  za,  is^ 
on  the  arc  OAj.  On  the  other  hand,  we  specify  the  quantities  E/b  and 
y  “ir  —  6 .  We  must  now  determine  whether  or  not  our  assumption  regarding 
the  location  of  za  can  restrict  possible  choices  for  the  parameters  E/b  and 
6.  In  particular,  we  recall  that,  as  is  true  for  the  quantity  we  must 
also  require  that  0  <  5a  <  5,  as  noted  in  Eq.  (45).  The  limiting  condition, 
corresponding  to  the  coincidence  of  the  apex  and  cavity  separation  point  at 
the  nose  of  the  profile  corresponds  to  £a  ■  0.  In  this  case,  the  smallest 
value  of  tj  for  any  choice  of  E/b  will  be  found  when 

aLn  b(1  •  C0S  6) 

On  the  other  hand,  by  inspection  of  Eqs.  (47)  and  (48),  we  see  that  the 
largest  values  of  f  and  g  are  found  for  ta  =  t^n.  The  values  of  F  and  G 
are  also  obtained  from  the  smallest  value  of  t2  because  tt  is  calculated 
from  t£  with  ^2  ”  namely: 
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t  E- 

t  b(  1  +  cos  6) 

Let  us  compare  the  values  of  F  with  f  and  G  with  g.  Suppose  that  E/b  is 
selected  so  that  cos  ta  »  cos  tt  ■  1,  sin  ta  •  ta  and  sin  t^  *  tf  Let 
6  <<  1.  In  this  case  we  can  see  that 


For  example,  if  6  =•  .1,  we  would  estimate  F/f  -  (20)^  and  G/g  ~  (20)^.  These 
estimates  imply  that  both  E/b  and  6  are  significantly  smaller  than  unity. 

In  the  applications  contemplated,  E/b  will  probably  be  less  than  unity 
although  6  might  conceivably  approach  or  exceed  unity.  Therefore,  we  shall 
consider  the  ratio. 


1  +  cos  6 
1  -  cos  6 


n  >  1  , 


which  permits  us  to  consider  roughly  the  ratios  of  F  to  f  and  of  G  to  g 
for  various  values  of  6.  In  particular,  we  can  solve  for  cos  6  and  obtain 


cos  6 


n  -  1 
n  +  1  * 


which  permits  us  to  plot  a  curve  of  n  *  (t  /t  )  versus  6  as  shown  in 

'min 

Fig.  8.  This  curve  illustrates  the  effect  that  the  choice  of  stagnation 


point  location  has  on  the  ratio,  n.  The  value  of  n  in  turn  gives  a  rough 
indication  of  how  large  the  ratios  F/f  and  G/g  will  be. 
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It  appears  for  most  cases  that  these  ratios  will  be  very  large  and 
one  need  not  solve  Eqs.  (49)  and  (50)  by  iteration.  Instead  one  can 
obtain  an  accurate  value  of  a  from 


tan  a  i  | 

and  he  can  then  determine  U/E^  from 

U  _  F 
g2  cos  a 


(49a 


(51a 


Should  cases  arise  in  which  Eqs.  (49a)  and  (51a)  are  not  accurate,  they  can 

be  used  advantageously  to  start  the  iteration.  In  order  to  illustrate  these 

points  and  in  order  to  show  a  profile  shape  derived  from  the  complementary 

solution,  we  have  prepared  the  following  numerical  example. 

We  started  the  calculation  by  selecting  6  *  70°  and  E/b  *  .01. 

Figure  8  shows  that  n  4  2.  The  values  t  I  *  .01520  and  t^  *  .00745 

a’min 

follow  from  the  formulae  for  these  quantities.  From  Eqs.  (47)  we  find  that 
F  *  18,004  and  G  ■  266.83.  Equations  (49a)  and  (51a)  lead  to  tan  o  *  .01482 
and  U/E^  »  18,006.  From  the  formulae  just  after  Eq.  (36)  we  have  for  the 
cavity  drag  due  to  a  profile  of  unit  chord 

E2 

CD  -  2ir  -g—  -  .00035 
c 


as  the  contribution  for  this  point -drag  profile.  The  value  of  ta  can  now 
be  found  from  Eq.  (50).  It  is  ta  -  1.556.  Equation  (43)  can  now  be  used 
to  find  that  £a  ■  69.61°.  Note  that  for  this  case  the  apex  is  almost 
coinc  ient  with  the  stagnation  point.  Because  Eq.  (35a)  shows  the 
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complementary  function  to  be  at  the  stagnation  point  this  result  is 
expected.  The  fact  that  the  apex  is  not  exactly  at  5  *  70°  is  due  only 
to  the  inclination,  a  *  .85°,  between  the  chord  line  and  the  x  axis. 

Once  U/e2  has  been  found,  the  values  of  xt  and  yt  can  be  found  from 
Eqs.  (47).  The  values  are  xt  =  .99989  and  yt  *  -  .01482.  Now  one  can 
use  the  conditions  that  in  Eqs.  (1),  (*t»yt)  **  (o=l,n*0)  and 
(xa.ya)  (o=0,n=0).  These  lead  to  a  system  of  equations  from  which 
oQ  and  nQ  can  be  eliminated  and  one  finds  that 

x  *  x  cos  a  =  x  -  sin  t  , 
at  t  a 

y  »  y  +  sin  a  =  y  +  cos  t  . 
a  t  a 

When  the  above  values  of  xt  and  o  were  used  in  these  equations,  the  values  of 
xa  and  ya  were  found  to  be  zero  to  within  five  decimal  places.  This  result 
is  consistent  with  the  location  of  £a  noted  previously.  Continuing  with 
Eqs.  (1),  we  can  use  the  fact  that  xa  =  ya  *  0  to  see  that  it  must  also 
follow  that  oQ  *  no  =  0.  Accordingly,  the  form  of  Eqs.  (1)  for  the  present 
calculations  is 

o  -  .99989x  -  .0 14796y  , 
n  =  .0 14796x  +  . 999989y  . 

The  next  phase  of  the  calculations  is  the  evaluation  of  the  equations 
for  the  wetted  surface  and  cavity  contours  in  accordance  with  Eqs.  (38), 

(39)  and  (40).  The  result  of  these  calculations  is  shown  in  Fig.  9.  In 
this  figure,  the  chordline  distance,  a,  has  been  labeled  as  X  and  the 
ordinate,  n,  has  been  labeled  as  Y.  Note  that  the  Y-scale  is  magnified 
five  times  compared  to  that  of  the  X-scale.  The  trailing  edge  of  the 
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M 


A 

I 


wetted  surface  is  at  X  *  1.  The  upper  surface  separation  point  is  at 
a  *  X  *  .240.  The  cavity  thickness  at  X  =  1  is  Y  =  T  3  .02980.  This 
point  is  marked  to  the  same  scale  as  the  X-scale  by  the  dot  and  the  line 
at  X  =■  1  in  order  to  give  an  idea  of  the  actual  thickness  of  this  example 
of  a  point -drag  profile.  Finally  we  can  calculate  the  value  of  Cc(l) 
in  this  example  in  order  to  compare  it  with  the  values  found  previously  for 
the  flat-plate.  It  is  found  from  the  value  of  t^  38  (E/b )/ [  (^c  +  1/ C c )/ 2  cos  y] 
at  a  3  1.  The  calculations  indicated  give  5^(1)  =  «3290  which  is  roughly  two 
times  the  value  from  Table  1. 

The  Eigensolution 

Our  desire  to  retain  as  much  simplicity  as  possible  in  the  preceding 
analysis  of  the  complementary  function  has  caused  us  to  place  the  point-drag 
singularity  at  the  stagnation  point,  and  that  is  why  we  call  it  the  simple 
eigensolution.  This  restriction  on  the  location  of  the  point  of  application 
has  allowed  us  to  show  that  such  a  solution  exists,  that  it  definitely  leads 
to  a  smoothly  rounded  profile  nose  and  that  it  will  cause  an  incremental 
thickening  of  the  cavity  depending  on  its  strength,  E.  Of  course,  we  need 
not  restrict  ourselves  to  the  stagnation  point  as  being  the  location  of  the 
point-drag  singularity. 

For  example,  suppose  we  choose  some  other  point  C  on  the  wetted  surface. 
Such  a  point  is  illustrated  in  Fig.  5  and  it  happens  to  be  located  between 
the  upper  cavity  separation  point  and  the  stagnation  point,  although  C  could 
just  as  well  be  at  some  other  wetted-surf ace  location.  The  main  idea  is  that 
now  0  *  Yc  at  the  location  of  the  point-drag  singularity  and  if  we  simply 
replace  Eq.  (34)  by  the  modified  expression, 


rm 


■i-V-V-. 


,  .  . . . 
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•c<«> 


b[cos  Yc  -  J  (c  +  £)] 


_ -  2Eg _ 

if  -iy 

b(5  -e  C)(?  “  e  c] 


b  sin  y 


C  -  e 


?  -  e 


-iy  J  * 
c 


(34a) 


we  still  have  a  function  which  satisfies  those  conditions  needed  for  a 

complementary  solution.  It  is  clear  that  in  the  £ -plane,  id  (0)  «  0  is  in 

c 

agreement  with  Eq.  (16).  Moreover,  when  5  is  real  cd  is  real  and  on  the 

c 

unit  circle  t  “0  everywhere  except  possibly  at  the  simple  poles, 

±iy 

C  “  e  C.  From  Eqs.  (22)  and  (23)  it  follows  that 


CD  a  2*k  v 

c 

as  before.  On  the  other  hand,  because  of  the  displacement  of  the  point  C 
away  from  0,  a  lift  force  is  produced  and  we  find  that  now 


81  bE  ,  1  +  '’c  . 

'L  “  “ifc-  *la  — 2 -  Sl” 

C 


Y  “  Y 


'c  _  8xbE  c 

2 - uT"  sin  — 


6+6  6-6 


,  (53) 


where  y  *  ir  -  6  and  yc  *  ir  -  6C  in  accordance  with  previous  convention.  The 
profile  chord,  c,  should  be  set  at  unity  in  Eqs.  (52)  and  (53).  Equation  (53) 
shows  that  Cl  *  0  when  6C  ■  6.  But  6C  <  6  when  the  point  C  moves  toward 
the  point  Aj  and  a  negative  lift  results.  In  the  limit  as  5C  +  0,  we  have 

bE  .26 
-  sin  Y  • 


-38- 


10  June  1985 
BRPrlhz 


If  C  is  between  0  and  A2  a  positive  lift  is  produced  and  in  the  limit  when  C 
is  coincident  with  A2  we  have 


8rbE 

Uc 


2  6 

cos  y  . 


If  one  were  to  let  E  be  negative  the  sign  of  the  foregoing  trends  with 
respect  to  C  would  be  reversed.  We  must  insist  however,  that  E  >  0  because 

Li 

this  function  produces  a  thickening  of  the  cavity  and  because  then 

u>'(0)  **  -  2E/b.  Similarly,  we  have  also  found  that  u»'(0)  -  -  2  sin  y 
c  o 

Thus,  the  effect  of  adding  uq  and  u>c  increases  the  net  drag.  Neither  of  these 
functions  can  act  to  reduce  it.  Accordingly,  we  shall  take  Eq.  (34a)  as  the 
appropriate  form  of  the  eigensolution  which  has  been  sought.  Since  both  of 
Eqs.  (34)  have  simple  poles  on  the  contour  | C |  ■  1,  they  can  be  thought  of  as 
elementary  solutions.  But  as  we  have  seen  in  the  case  of  the  simple  eigen¬ 
solution,  the  pole  at  8  *  y  does  not  lead  to  an  unbounded  value  of  q  /U. 

c  c 

Indeed,  0  <  q  /U  <  1  in  the  neighborhood  of  (5  *  ±  y  ,  which  can  be  taken  to 
c  c 

be  at  any  point  on  the  wetted  surface  when  8  *  +  y  . 

c 


Some  Profile  Geometry  and  the  Flow 

In  any  inverse  design  procedure,  one  starts  the  calculation  by 
prescribing  the  pressure  distribution  or  the  magnitude  of  the  velocity 
along  the  periphery  of  the  profile.  Of  course,  that  is  almost  like 
having  the  solution  at  the  outset;  but  not  quite,  because  one  does  not 
know  the  relationship  between  points  along  the  hydrofoil  arc  length,  s, 
and  corresponding  points  on  the  semi -circle  in  the  5  plane. 
Nevertheless,  referring  to  Fig.  7,  we  can  measure  the  arc  length  s  from 


the  point  A^  at  the  cavity-trailing  edge  separation  point.  Then  the 
arc  length  Increases  from  A.^  until  one  reaches  the  stagnation  point  0  in 
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in  Fig.  7.  At  0  the  arc  length  will  be  designated  by  sq.  Continuing  along 

the  periphery,  one  rounds  the  nose  of  the  foil  and  arrives  at  the  separation 

point  A^.  At  this  point,  the  arc  length  is  s^.  Finally  we  proceed  along 

the  upper  surface  of  the  cavity  until  we  arrive  at  the  point  s2»  a  distance 

T  directly  above  the  trailing  edge,  A2*  The  distance  T  is  measured 

perpendicularly  to  the  profile  chordline.  Thus,  it  is  parallel  to  the  n 

axis  in  Fig.  7.  Clearly,  0<s  <s,<s„. 

o  1  2 

A  schematic  diagram  showing  the  flow  speed  on  the  wetted-surf ace  arc 

and  the  upper  surface  of  the  cavity  is  illustrated  in  Fig.  10.  This  figure 

shows  rather  clearly  that  the  designer  does  not  have  as  much  freedom  with 

regard  to  the  pressure  distribution  prescription  as  he  might  wish.  For 

example,  we  know  that  |q/u|  «  1  at  s  ■  0,  s  «  Sj  and  in  the  interval, 

S1  <  S  <  S2’  Moreover  we  take  the  values  of  q/U  <  0  in  the  interval 

0  <  s  <  sq  because  the  flow  direction  on  the  wetted  surface  points  in  the 

direction  opposite  to  the  positive  sense  of  s.  Between  s  <  s  <  s_  the 

o  2 

opposite  situation  holds  and  we  count  q/U  >  0  in  this  interval.  At  no 
point  in  the  flow  can  q/U  >  1. 

Since  we  prescribe  the  magnitude  of  q/U  everywhere  this  is  the  same  as 
prescribing  the  lift  coefficient  C  .  The  prescribed  value  of  C  can  be 
used  to  fix  the  relative  position  of  the  stagnation  point  Sq  with  respect 
to  point  s2  if  the  distribution  q(s)/U  is  not  too  firmly  fixed.  One  need 
only  use  the  well-known  Kutta-Joukowsky  formula. 
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where  c  is  the  profile  chord  and  T  is  the  circulation.  But  we  can  write 


where 


r  -  (a2  -  ax)u 


Ms0)  -/  V~ds 


4(so»s2)  “  /  +  («2  "  sl)  * 


Generally  one  will  need  to  solve  for  the  ratios  s  /s„  and  s,/s,  numerically, 
starting  the  iteration  by  supposing  that  possibly  A^  and  that  part  of  A^ 
between  sq  and  s^  are  triangular  areas.  Indeed,  we  will  see  below  that  the 
area  A^  must  be  rather  closely  specified  in  advance  so  that  the  chief  freedom 
to  be  exercised  by  the  designer  is  associated  with  A^ .  The  thought  behind 
these  observations  is  that  whatever  the  approach,  the  design  CL  and  the 
distribution  q(s)/U  must  be  consistently  prescribed. 

We  now  turn  to  the  properties  of  the  flow  in  the  neighborhoods  of  s 


and  Sj  because  at  sq  the  flow  on  the  wetted  surface  is  a  stagnation  flow  and 

at  s^  it  is  constrained  by  the  requirement  for  smooth  separation.  For 

simplicity's  sake  we  shall  shift  our  reference  point  from  to  the  point  0 

and  measure  the  arc  length  s  from  0.  This  means  that  the  distance  along  the 

arc  from  0  to  A,  is  s,  -  s  and  we  will  normalize  all  intermediate  distances 
1  1  o 


s  by  writing  x  ■  s/fs,  -  s  ), 

1  o 
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Then  it  is  known  for  a  potential -flow  stagnation  point  on  a  flat  wall 
that  the  streamlines  are  equilateral  hyperbolas  having  their  separatrices 
as  the  straight  wall  and  the  normal  to  the  wall  through  the  stagnation  point. 
It  is  also  known  from  linearized  theory  [9] ,  and  it  can  also  be  shown  for  the 
exact  theory,  that  if  x  -  1  is  a  separation  point  then  C  (x)  ~  /  1  -  x.  A 

p 

specific  example  of  this  general  behavior  can  be  seen  by  referring  to  the 

various  formulae  for  the  flat-plate  velocity  and  pressure  distributions 

given  above  in  the  discussion  surrounding  u>  .  In  particular,  it  is  easily 

o 

seen  that  dz/dij  *  0  at  the  separation  points  5  *  1.  Consequently  curves 

of  q/U  and  will  have  vertical  tangents  at  the  separation  points.  In 
response  to  these  requirements  we  shall  consider  a  one-parameter  family  of 
speed  distributions  on  the  forward  part  of  the  wetted  surface;  namely, 

q(x)/U  *  x  [l  -  a  Vi  -*2  (1 -V7  -  X2)]  ,  (0  <  X  <  1)  .  (54) 

Although  this  is  a  rather  special  class  of  velocity  distributions,  it 
probably  contains  as  much  generality  as  one  usually  needs  because  most  often 
the  arc  length  from  sq  to  s^  is  very  small  compared  to  the  total  arc  length 
s^.  Therefore  the  curves  of  q(x)/U  will  be  very  steep  and  the  shape  of  one 
choice  of  distribution  would  hardly  be  discernible  from  some  other  choice, 
provided  that  the  required  conditions  at  x  ■  0  and  x  *  1  are  satisfied. 

Plots  of  Eq.  (54)  are  given  for  a  highly  stretched  length  scale  in  Fig.  11. 

This  specification  and  cavity  surface  velocity  for  s^  <  s  <  s^  makes 
the  entire  pressure  distribution  on  the  upper  surface  of  the  profile  fairly 
well  defined.  In  order  that  these  results  can  be  used  conveniently  in  the 
process  of  reconciling  the  design  C  with  q(s)/U,  Eq.  (54)  should  be 

Li 

expressed  in  terms  of  the  arc-length  coordinates  defined  in  Fig.  10. 
Normalizing  all  distances  by  s ^  we  can  write 
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8,  8  “  8 

x  -  ( - 1 - )( - 2.)  . 

8l-8on  8  2 


Having  prescribed  the  velocity  distribution  for  points  on  the  periphery 
of  the  profile,  one  can  now  use  Eq.  (20)  in  order  to  derive  the  circle-to- 
profile  correlation  between  s^  and  s^  albeit  with  some  degree  of  ambiguity 
until  the  entire  problem  has  been  formulated  and  solved.  For  this 
correlation,  we  will  write  Eq.  (20)  in  the  form, 


q(s) 


2  —  (cos  y  -  cos  B)sin  8  d8  *  d[(cos  y  -  cos  B)2]  ,  (56) 


where  Eq.  (54)  defines  q(s)/U  in  the  interval,  s  <  s  <  s  ,  and  in  the 

o  1 

interval,  s^  <  s  <  s^,  we  have  q(s)/U  »  1.  Moreover,  at  this  point  we  know 
that  when  s  ■  s^,  8  *  Y  and  when  s  ■  s^,  0  *  if .  From  Eq.  (55)  we  can  put 


(s,  -  s  )dx  ■  ds 
^  1  oJ 


and  we  can  integrate  Eq.  (56)  from  sq(x  ■  0)  to  x  and  the  right-hand  side  of 
Eq.  (56)  from  8  -  Y  to  8.  Therefore  the  integral  of  Eq.  (56)  is 


(81  -s0)[—-a( 


2  2  ' 
x4  _r  1  -  (1  -  xZ) 


-  f-  +  f-)]  -  Jj-  (cos  Y  -  COS  b)2  . 


When  8  *  if  and  x  «  1,  Eq.  (57)  becomes 


(S1  ~  80)(l  -  f)  -  2  jj-  (1  +  cos  y)‘ 
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Proceeding  Co  Che  arc  <  s  <  along  Che  upper  surface  of  Che  cavicy 
one  may  sCarC  ac  Che  poinC  0  -  it  on  Che  circle  |c|  -  1  and  move  along  Che 
negacive  real  axis  Coward  £  ■  -  r  ■  0.  The  velociCy  q/U  ■  1  everywhere 
along  chis  arc  and  Eqs.  (19)  can  be  used  Co  find  Che  differencial  arc-lengch 
correlaCion, 

d8  *  -  jp  [j  (7  +  r)  +  cos  *](■£■  -  r)~  ,  (59) 

where  we  have  Caken  Che  posicive  square  rooC  in  view  of  Che  face  Chac  r  has 
been  defined  as  a  posicive  quanCiCy.  SCarCing  Che  inCegraCion  ac  s  -  s^  and 
r  *  1,  we  find 

2 

8  "  Sj  *  ]j“  {(7  “  l)[^  (7  +  0  +  cos  y]  -  (1  -  r)[-£  (1  +  r)  +  cos  y]}  .  (60) 

If  r  «  r^  when  s  ■  s^  we  can  sum  Chis  special  case  of  Eq.  (60)  and  Eq.  (58) 
in  order  Co  geC  expressions  for  in  Che  profile  and  c  planes.  As  implied 
from  Che  ouCseC,  Che  consideraCions  of  Chis  secCion  apply  Co  Che  inverse 
problem  as  a  whole  and  we  could  exCend  the  above  consideraCions  to  deduce 
a  correlation  like  (57)  for  the  interval  between  k^  and  0.  We  shall  not 
complete  Che  analysis  here,  however. 

The  Role  of  the  Eigensolution 

Returning  to  the  eigensolution  as  given  by  Eq.  (34a),  we  recall  that 

±iYc 

except  for  the  simple  poles  at  ;  ■  e  ,  the  imaginary  part  of  u>c  is  zero 
everywhere  on  the  unit  circle  and  the  real  part  suffers  a  jump  of  magnitude 
v  at  t  on  the  unit  circle.  On  the  other  hand,  we  have  found  that  the 
complex  velocity  has  an  isolated  essential  singularity  at  these  poles  (see 


p.  22  above).  Moreover,  because  of  this  fact  we  observed  that  q  /U  is 

c 

bounded,  being  in  the  interval  [0,1].  Therefore  by  appealing  to  a  theorem 

of  Weirestrass*,  we  shall  assign  the  value  q  /U  ■  1  at  the  points  y  ■  ±  1 

c  c 

on  | c|  ■  1.  We  have  also  seen  from  an  exploration  of  the  flow  due  to  an 

isolated  eigensolution  that  the  resulting  profile  is  smooth  and  that  its 

wetted  surface  is  smooth.  When  the  flat  plate  function  uiq  and  eigenfunction 

are  combined,  the  resulting  profile  will  also  have  a  stagnation  point  at 

6  ■  y,  as  we  shall  discuss  below. 

As  we  have  remarked  above,  our  plan  is  to  write  the  logarithmic 

hodograph  as  the  sum  of  the  flat  plate  function  u  ,  the  eigenfunction  u 

o  c 

and  a  regular  function  Thus 


qrt  q„  q, 

o  c  1  - 


(0(C)  -  o.o  +  (0c  +  *  0o  +  8c  +  eL  +  iin  [—  —  — ]  . 


(61) 


Consequently, 


1  -  qoqcql 


U 


(62) 


U 


The  composite  function,  q/U,  is  prescribed  from  the  outset.  Equation  (54) 
which  has  been  defined  over  the  arc  length  s  in  the  interval  [s  ,8^] 
corresponding  to  the  normalized  variable  x  in  the  interval  [0,1],  provides 


*See  for  example,  Copson,  E.  T.,  Theory  of  Functions  of  a  Complex  Variable 
Oxford  University  Press,  1946,  p.  81;  or  Tichmarsh,  The  Theory  of  Functions 
2n<*  Edition,  Oxford,  1949,  pp.  93-94.  ~  ~ 


f*t*V 


an  example  of  such  a  composite  function.  The  flat  plate  function,  q  /U,  has 

o 

been  worked  out  on  p.  17.  Here  it  will  be  written  as 

%  _ 2(-  cos  3  +  cos  y) _ 

U  2  2  *  (63) 

(-  cos  8  +  cos  y)  +  (sin  3  +  sin  y) 

where  the  absolute  sign  in  the  numerator  from  p.  16  has  been  replaced  by 
ordinary  brackets  and  we  have  made  other  changes  which  apply  in  the  interval 
Y  <  8  <  »•  Since  q^/U  *  1  on  |c|  ”1  the  function  q^/U  follows  from  Eq.  (62) 
rewritten  as  q^/U  *  (q/U)/ (qQ/U) .  For  the  example  at  hand,  the  result  of  this 
transposition  can  not  be  used  until  we  transform  Eq.  (63)  from  the  c  plane 
into  the  arc  length  in  the  interval  [s^.s^]  with  the  help  of  Eq.  (57).  In 


particular,  we  find  that 


2  2  2 
sin  8  »  1  -  (R  -  cos  y)  *  sin  Y  +  2R  cos  y 


-  R2  , 


where 


„2  U  t  ux 

E  *TT  (si  ~s„% 

D 


(1  -  X^) 

3 


2, 3/2  2  4 


Equation  (64)  enables  us  to  write  Eq.  (63)  in  terms  of  arc  length 


R  +  sin  Y  +  V  sin  Y  +  2R  cos  y 
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Then  it  follows  from  Eq.  (62)  that 

if*  *  {(sl  "  s0)x[1  "  1  "  x2  -  V  1  -  x2)]}  *■ 

{2R/[R2  +  sin  y  +  V  sin2  y  +  2R  cos  ?  -  R2]}  .  (66) 

As  is  true  of  q  /U  and  q  / U,  the  function  q,/U  when  represented  in  the  £  plane 
o  c  1 

can  be  continued  to  the  arc  of  the  unit  circle  in  the  lower  half  of  the 

C -plane  in  accordance  with  the  formulas  at  the  bottom  of  p.  13. 

Now  that  the  well  behaved  function  q^/U  has  been  found  along  the  arc  OA^, 

we  may  return  to  Eq.  (61)  and  note  that  the  real  part  of  co(c)  is  not  yet 

known  completely  even  in  this  interval.  It  will  be  known  however,  if  we  can 

find  9j.  But  its  complex  conjugate  x^  has  just  been  found.  Of  course,  we  do 
not  know  x^  at  all  points  of  the  unit  circle  because  in  this  example  we  have 
have  not  prescribed  the  entire  pressure  distribution  on  the  wetted  surface, 
although  the  way  in  which  this  can  be  done  is  certainly  clear.  Once  that 
step  has  been  carried  out,  x^  will  be  known  on  the  unit  circle  and  its 
conjugate  9^  can  be  calculated  using  the  customary  representation  of 
u)^(0  as  a  Laurent  expansion,  as  employed  by  Yoshihara  [8]  for  example. 
Procedures  which  are  preferable  for  practical  engineering  calculations  are 
usually  based  on  the  Poisson  integral  formula  or  related  methods.  Examples 
of  interest  for  the  present  inverse  theory  are  given  by  Theodorsen  and 
Garrick  [16]  and  by  Parkin  and  Peebles  [17],  among  others.  Further 
discussion  of  these  matters  is  beyond  the  scope  of  this  report. 


-v-:, *>. £-.v vS-Sv-i 
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In  the  course  of  designing  a  profile,  the  quantities  which  one  specifies 
from  the  start  are  the  design  lift  coefficient  C^,  the  pressure  (velocity) 
distribution  on  the  wetted  surface,  the  cavity  thickness  at  the  trailing  edge 
and  perhaps  the  separation  point  of  the  upper  surface  of  the  cavity  near  the 
nose.  The  purpose  of  the  eigensolution  is  to  provide  the  necessary  degrees 
of  freedom  which  will  permit  the  control  of  the  cavity  geometry  as  indicated. 
Added  degrees  of  freedom  can  be  incorporated  in  the  prescribed  pressure 
distribution  if  a  parametric  approach  is  used.  It  might  be  possible  to  lower 
the  cavity  drag  somewhat  by  adjusting  them  although  the  prescription  of 
cavity  thickness  is  probably  more  important  in  this  regard.  The  outcome  of 
the  design  process  will  be  coordinates  of  the  profile  and  cavity  shape, 
including  the  separation  point;  the  attack  angle,  a;  and  the  drag  coefficient 
Cp.  The  eigensolution  strength,  E/b,  will  also  be  determined  in  the  course 
of  the  calculations. 

Conclusions 

The  chief  finding  of  this  paper  is  that  one  can  construct  many  singular 
eigensolutions  for  the  exact  inverse  problem  of  two-dimensional  cavity  flow  at 
zero  cavitation  number.  From  among  these,  we  have  chosen  that  single 
eigensolution  which  provides  the  correct  branching  of  the  flow  at  its 
singularity,  it  also  appears  to  offer  the  least  disturbance  to  the  upstream 
flow  field  inclination  of  any  cavity  flow  which  does  not  already  include  a 
point-drag  solution  as  one  of  its  elements.  This  particular  choice  also  seems 
to  offer  the  greatest  analytical  convenience.  The  physical  conditions 
satisfied  by  this  eigensolution  are: 
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(1)  At  points  on  the  cavity  and  on  the  wetted  surface  of  the 
profile,  the  flow  velocity  is  equal  in  magnitude  to  the 
free-stream  velocity.  Consequently,  except  for  the 
singular  point  the  pressure  coefficient  is  zero  on  the 
wetted  and  cavity  surfaces. 

(2)  The  point-drag  solution  vanishes  at  infinity,  but  it 
does  have  a  bounded  essential  singularity  on  the  wetted 
surface  and  it  produces  no  singular  velocities  or 
pressures  in  the  flow. 

(3)  This  function  produces  no  additional  flow  inclination 
on  the  entire  upstream  stagnation  streamline. 

A  specific  example  of  the  flow  geometry  represented  by  an  isolated 
eigensolution  has  been  given  above  to  show  how  this  function  can  produce 
round-nosed  profiles.  In  general,  it  is  found  that  the  point-drag  solution 
produces  a  widening  of  the  cavity  which  is  directly  proportional  to  its 
strength.  An  incremental  cavity  drag  accompanies  this  widening  and  this 
drag  is  proportional  to  the  square  of  the  eigensolution  strength.  No  lift 
is  produced  by  the  point -drag  function  when  its  location  coincides  with 
that  of  the  stagnation  point  on  the  profile  surface.  In  contrast  to  the 
linearized  theory,  the  complementary  function  singularity  need  not  be  at 
the  stagnation  point.  In  these  cases,  the  incremental  cavity  drag  is 
not  changed  from  its  value  when  the  singularity  is  at  the  stagnation  point. 
But  when  the  singularity  is  located  between  the  stagnation  point  and  the 
upper  separation  point  a  negative  incremental  lift  is  produced.  If  the 
singularity  is  on  the  lower  surface,  downstream  of  the  stagnation  point. 
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a  positive  lift  increment  is  found.  Whether  or  not  it  is  better  to  position 
the  eigensolution  at  some  point  on  the  profile  instead  of  at  the  stagnation 
point  remains  to  be  studied. 

As  a  result  of  these  findings,  it  appears  that  an  eigensolution  exists 
for  the  nonlinear  theory  of  cavity  flow  at  zero  cavitation  number  and  that 
it  is  now  most  likely  that  a  similar  eigensolution  can  be  found  for  such 
fully  cavitating  flows  at  cavitation  numbers  which  are  greater  than  zero. 

The  results  found  so  far  suggest  that  the  nonlinearized  theory  and  the 
linearized  theories  parallel  one  another  very  closely  as  far  as  the  nature 
of  their  point  drag  solutions  are  concerned.  But  the  present  results  exhibit 
some  featues  which  are  lost  in  the  process  of  linearization. 
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Figure  1 


Profile  and  cavity  geometry  in  the  physical 
or  z-plane  at  zero  cavitation  number. 
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Figure  2.  Profile  and  cavity  surfaces  at  zero  cavitation  number 
in  the  plane  of  the  complex  potential,  F  *  (f>  +  iijj. 
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Figure  3.  The  mapping  W  =  v^F  maps  the  flow  outside  the  cut  in 
the  F  plane  into  the  upper  half  W  plane. 
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Figure  4.  The  wetted  surface  lies  between  -1  and  +1  in  the 
complex  Z  plane. 
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Figure  5.  The  Joukowskl  transformation  maps  the  flow  from  the 
upper  half  of  the  Z-plane  Into  the  interior  of  the 
upper  unit  semi-circle  with  the  point  at  infinity 
at  the  origin. 
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vertical  and  horizontal  axes  have  the  same  scale. 


Figure  10.  Schematic  diagram  showing  a  prescribed  velocity  distribution 
on  a  profile  and  cavity  arc. 
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Figure  11.  A  one-parameter  family  of  velocity  and  pressure  distributions 
over  the  nose  arc  of  a  hydrofoil. 
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